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Abstract. The upcoming XMM Large Scale Structure Survey (XMM-LSS) will ultimately provide a unique map- 
ping of the distribution of X-ray sources in a contiguous 64 deg^ region. In particular, it will provide the 3- 
dimensional location of about 900 galaxy clusters out to a redshift of about 1. We study the prospects that this 
cluster catalogue offers for measuring cosmological parameters. We use the Press-Schechter formalism to predict 
the counts of clusters and their X-ray properties in several CDM models. We compute the detection efficiency 
of clusters, using realistic simulations of XMM X-ray images, and study how it differs from a conventional flux 
limit. We compute the expected correlation function of clusters using the extended halo model, and show that 
it is expected to evolve very little out to z ~ 2, once the selection function of the survey is taken into account. 
The shape and the amplitude of the correlation function can be used to brake degeneracies present when cluster 
counts alone are considered. Ignoring systematic uncertainties, the combination of cluster counts evolution and 
of the correlation function yields measurements of Qm, o"8 and F with a precision of about 15%, 10% and 35%, 
respectively, in a ACDM model. This combination will also provide a consistency check for the ACDM model, 
and a discrimination between this model and the OCDM model. The XMM-LSS will therefore provide important 
constraints on cosmological parameters, complementing that from other methods such as the Cosmic Microwave 
Background. We discuss how these constraints are affected by instrumental systematics and by the uncertainties 
in the scaling relations of clusters. 

Key words. X-rays: galaxies: clusters; Galaxies: clusters: general; Cosmology: cosmological parameters; Cosmology: 
large-scale structure of Universe; Surveys 



1. Introduction 

Clusters of galaxies are the most massive bound objects in 
the Universe and provide a pow erful cosmological probe 
(see e.g. Borgani & Guzzo 2001 for a review). In particu- 



lar, the number counts of clusters and its evolution yield a 
robust measure of both the amplitude of the matter power 
spectrum and of the geometry of the universe (e.g. Oukbir 
Eke et al 



& Blanchard 1997 



1998; Viana & Liddle 1999 



Kitiyama & Suto 1997). The spatial correlation function 
of clusters quantifies the clustering of these objects and 
yields complementary constraints on cosmology (e.g. Mo, 
Jing & White |l996t Suto et al. ^OOOj R obinson pOOC ; 



Moscardini et al. 2000; Collins et al. 2000) 



In this paper, we explore the prospects of measuring 
cosmological parameters with the upcoming XMM Large 
Scale Structure Survey (XMM-LSS; Pierre |2000| ). This 



survey consists of 10 ksec exposures of an 8 x 8 deg re- 
gion with the XMM-Newton observatory, along with an 
extensive follow-up programme in the optical, IR and ra- 
dio bands. In particular, it will provide the 3-dimensional 



location of about 900 clusters out to a redshift of about 
1. Thanks to its uniform sensitivity across a contiguous 
region, this survey thus provides a unique database to 
measure the evolution of both the number counts and the 
correlation function of clusters. 

To study how the clusters found in XMM-LSS can con- 
strain cosmological models, we use the Press & Schechter 
( 1974 ) formalism to predict the expected cluster counts 
in the survey. This is done using the selection function of 
the survey derived from detailed simulations of cluster de- 
tections in XMM-Newton images (see Valtchanov, Pierre 
& Gastaud pOOl| , VPG). Using the Mo & White ( [l996| ) 
formalism, we compute the expected spatial correlation 
function for the detected clusters. We then study how the 
cluster counts and correlation function, taken together, 
constrain cosmological parameters. 

Ou r analysis extends the work of Moscardini et al. 
( ^OOOl) who considered the expected cluster counts and 
correlation function for XMM-LSS. They however as- 
sumed a simple flux limit, rather than the more real- 
istic selection function which we consider. In addition. 
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they did not compute the cosmological constraints from 
a joint measurement of the cluster counts and correla- 
tion function with XMM-LSS. Our results also comple- 
ment the analysis of Romer et al. ( ^001 ) who studied the 
cosmological dependence of cluster counts in the planned 
Serendipitous XMM Cluster Survey. They are also related 



to the work of Haiman al. (2001) and Holder et al. (2001) 
who studied the cosmological constraints which can be de- 
rived from the evolution of cluster counts with a dedicated 
wide-angle X-ray mission. 

The paper is organized as follows. In §^ we summa- 
rize the characteristics of the XMM-LSS. In we de- 
scribe the simulations for cluster detection and derive the 
cluster selection function. In §^ we compute the expected 
cluster counts using the selection function combined with 
the Press-Schechter formalism. In ^ we compute the cor- 
relation function for this cluster sample and show how it 
constrains cosmological parameters. The effects of system- 
atic uncertainties on these constraints are discussed in 
Our conclusions are summarized in §0. 



2. The XMM Large Scale Structure Survey 

The XMM-LSS Survey is a unique medium-deep cluster 
survey combining X-ray observations with an extensive 
optical, IR, and radio follow-up programme (Pierre 2000| ). 
The survey geometry - coverage and depth - was chosen to 
allow the measurement of the cluster two-point correlation 
function, with better than 15% error on the correlation 
length, in two redshift intervals between z = and 1. 

The position of the 8x8 deg^ surveyed area on the 
sky (a = 2'*20™, S — —5°) is at a sufficiently high galactic 
latitude (« —60°) in a region of moderate galactic ab- 
sorption and without known bright X-ray sources. It will 
be covered by 24 x 24 partially overlapping XMM point- 
ings with individual exposure times of 10 ks, reaching a 
sensitivity of 3 x 10^^^ erg s^^ cm^^ in the [0.5 — 2] keV 
band for point sources, or of about 5 x 10^^'^ erg s^^ cm^^ 



for cluster-like extended sources. Down to this limit, some 
300 objects (mainly QSOs) per deg^ are expected accord- 
ing to the latest deep surveys (e.g. Hasinger et al. 2001 



Giacconi et al. 2001), and a total of about 10-15 clusters 



per deg out to 2: ~ 1. The survey is also well suited to 
probe the existence of massive clusters within the impor- 
tant 1 < z < 2 redshift range. At the time of writing, 6 
deg^ with XMM have been allocated in guaranteed time 
and guest observer time. The rest of the survey is subject 
to an ongoing application and reviewing process. 

An extensive multi-wavelength follow-up programme 
has been undertaken by the XMM-LSS consortiurrj^. 
Special care is given to the optical identification of the X- 
ray sources: deep multi-color imaging of the entire region 
will be performed by the Canada-France-Hawaii Legacy 



Survey^ and subsequent redshift measurements by the 
VIRMOS/VLT instrument and other large telescopes to 
which the consortium has access. The main priorities are: 
(1) identification and redshift measurement of all X-ray 
clusters between < z < 1, (2) NIR observations of dis- 
tant (z > 1) cluster candidates and, subsequently, deter- 
mination of their redshift, (3) serendipitous spectroscopic 
observations of the X-ray QSOs, in order to study their 
clustering properties within the deep potential-well net- 
work traced by the clusters. 

3. Simulations 

In order to estimate the detection probability of clusters 
in the survey, we performed a series of X-ray image simu- 
lations. While a detailed description can be found in VPG, 
we first review here the main features of the simulations. 
We then show how they can be used to derive the selection 
function for the survey. 

3.1. Cluster Detection in Simulated Images 

The simulations reproduce the main characteristics of 
the XMM-EPIC instruments, such as the Point-Spread 
Function (assuming circular symmetry for the PSF shape; 
see VPG) and vignetting as a function of energy and 
off- axis angle. These were parametrized using the latest 
available on-flight calibrations. For the diffuse and parti- 
cle ba ckground, we have used the data from Watson et al. 
( ^00 ll) . 

The point-like sources were laid at random inside the 
field-of-view. Their fluxes were drawn from t he log N — 
logS data in Lockman Hole (Hasinger et al. 2001) and 
Chandra deep field south (Giacconi et al. 2001). Their 
spectrum was taken to be a power law with photon index 
2. 

Clusters of galaxies were modeled as spherically sym- 
metric objects assuming a /3-profile (e.g. Cavaliere & 
Fusco-Femiano 1976 ) with fixed core radius rc — 125 h^^ 
kpc and slope /3 = 0.75 (see §H for a discussion of 
the effect of a varying core radius on our predictions). 
A thermal plasma spectrum (Raymond & Smith 1977) 
was assumed. The spectrum was normalized using the 
non-evolving luminosity-temperature (L — T) relation of 
Arnaud & Evrard (1999). When generalized to arbitrary 
cosmological model, it is given by 

where L is the bolometric luminosity, T is the X-ray tem- 
perature, and Dl and I?L,EdS are the luminosity distances 
in the desired and Einstein-de Sitter cosmological models, 
respectively. 

We took the neutral hydrogen column density to be 
Nh = 5 X lO^*' cm~^ and element abundances Z = O.3Z0. 



Oflicial web page of the consortium: 



http:/ /vela.astro. ulg.ac.be/themes/spatial/xmm/LSS/ 



^ |http:/ /cdsweb.u-strasbg.fr:2001/Instruments/Imaging/ 
Megacam 
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Using XSPEC (Arnaud 1996), we calculated the total ex- 
pected count rates for the extended and point-like sources 
for the three XMM EPIC instruments^ with thin filters in 
[0.5-2] keV energy band, for an integration time of 10 ks. 

For a given temperature and redshift, we placed 13 
clusters with centers on a grid inside the inner 10' radius of 
the field-of-view. This number was chosen to maximally fill 
this region while avoiding overlaps. The grid pattern can 
be seen in the right-hand panel of Figure |l|. The detection 
was performed in the same way as in VPG, i.e. using multi- 
scale (wavelet) filtering assuming Poisson noise statistics 
(Starck fc Pie rre pJ9^ ), followed by SExtractor (Bertin & 



Arnouts 1996) detection and classification. As was pointed 
out in VPG, this is currently the most suitable method to 
detect, characterize and classify extended sources in XMM 
images. The raw photon and wavelet filtered images for 
clusters with T = 3 keV at a redshifts of z = 1 and 1.5 
are shown in Figs. |l]and| 

To cross-identify the detected objects with the input 
clusters we have used a searching radius of 12". If a cor- 
respondence is found, we perform a classification based 
on the half-light radius and the stellarity index, to de- 
termine whether the object is extended. (For the choice 
of the searching radius and the classification criteria, see 
VPG). A cluster is finally considered to be detected if the 
positional and the classification criteria are obeyed. This 
procedure is close to the planned analysis of the incom- 
ing XMM data, which will make use of the multi-color 
optical data to confirm the existence of the X-ray cluster 
candidates. 



3.2. Selection Function 

The selection function is a convenient way of quantifying 
the detection sensitivity of a survey. It was computed for a 
number of existing and future cluster surveys (eg. Romer 
et al. 2001 for the XMM serendipitous Cluster Survey; 



Henry et al. 2001 for the NEP ROSAT survey; Adami 



et al. 200C for the SHARC survey) using both analytical 
estimates and image simulations. 

Here, we estimate the selection function of the XMM- 
LSS, by performing a set of 10 simulations for a set 
of temperatures (T = 3,4,5,6,7,8,9 keV) and redshifts 
{z = 0.6, 1.0, 1.5, 1.8, 2.0). For each value of T and z, the 
selection (j){T, z) of clusters for the survey was then calcu- 
lated by comparing the number of detections (and correct 
classifications) N^ctiT, z) to the number of input clusters 
Nin{T, z), so that 



0(r,z) = iVdct(T,z)/iVi„(r,z). 



(2) 



Fig. ^ shows the resulting selection function, which gives 
the probability for a cluster with temperature T and red- 
shift z to be detected and classified as an extended object 
in the survey catalogue. Approximately, 90% of all clus- 
ters with T > 3 keV are detectable out to z ~ 0.6. The 



selection function is close to about 1 for T > 2 keV and 
z < 0.5 (not shown). Since we are only interested in clus- 
ters (and groups), we set the selection function to for 
T < 2 keV at all redshifts. This corresponds to a mini- 
mum bolometric luminosity of about 1.2 x 10^^ h~'^ erg 
s^^ (see Eq. Q) and can thus be easily implemented in 
practice. As can be seen on the figure, low temperature 
clusters become progressively harder to detect as the red- 
shift increases (compare Figs. ^ and For example, at 
z — 2 only clusters with T > 6 keV yield a completeness 
better than 90%. 

Our results are consistent with that of Romer et al. 



^ See e.g. the XMM-Newton User's Handbook: 
http: / / xmm.vilspa.esa.es / user /uhb / xmm_uhb.html 



( 2001 ; figure 5) who found that essentially all clusters with 
with T > 2 keV (T > 4 keV) and z <, 0.6 {z <, 1.0) can 
be detected in a survey with a sensitivity comparable to 
that of XMM-LSS. Note that our selection function corre- 
sponds to the average sensitivity over the inner 13 arcmin 
of a single XMM-EPIC pointing. The effects of sensitivity 
variations produced by vignetting and overlapping point- 
ings are discussed in ^ 

It is instructive to compare our selection function to 
that corresponding to a constant flux limit, as assumed 
in many previous studies. The flux 5'[o.5_2] of a cluster in 
the [0.5 — 2] keV band at redshift z can be derived from 
its temperature using the L — T relation (Eq. |^]), the 
Raymond-Smith spectrum and the luminosity-distance re- 
lation. Using this correspondence, we can express the se- 
lection function in terms of the flux rather than tem- 
perature. The resulting selection function ^(>S'[o.5_2] , 2) is 
shown in Fig. ^ Our selection function clearly does not 
exactly correspond to a sharp flux limit, but is instead 
a smooth function of the flux. For latter comparison, we 
note that a detection probability of 50% corresponds ap- 
proximately to a flux limit of 5'[o.5-2] = 10^^"* erg s~^ 
cm~^ at the redshifts of interest. 

In the following, we will use our selection function 
which includes most of the relevant instrumental and ob- 
servational limitations. For comparison, we will also con- 
sider the selection function corresponding to the above 
flux limit, along with that corresponding to temperature 
limit T > 2 keV. The latter selection is useful to study the 
effect of the removal of small clusters on our predictions. 



4. Cluster Counts 

We first compute the expected cluster number counts in 
the survey. This is done using the cluster selection function 
derived in the previous section combined with the Press- 
Schechter formalism. We first briefly review the main as- 
sumptions involved in our calculation of the mass function 
and of the temperature function of clusters. We then com- 
pute the expected projected number of clusters on the sky, 
as function of redshift. Finally, we show how the resulting 
redshift distribution constrains cosmological parameters 
within CDM models. 
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Fig. 1. Simulated extragalactic XMM-LSS raw image (left) and the corresponding wavelet filtered image with Poisson 
noise model and 10~^ (~ 4ct) significance (Starck & Pierre 1998). The energy band is [0.5 — 2] keV, the exposure 
time lOks, and the three XMM instruments are added together (MOSl, M0S2 and pn). Point-like sources follow the 
observed log N — log S relation and the extended sources (clearly visible on the wavelet filtered image) are clusters of 
galaxies with Tx = 3 keV at z = 1 (see text for details). 
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Fig. 2. Same as Fig. ^ but with clusters at z = 1.5. 



4.1. Mass Function 

The Press-Schechter formalism provides an analytic ex- 
pression for the a bundance of dark matter halos (Press 
& Schechter 1974 ). At a given redshift z, the differential 
number of dark matter halos of mass M per unit comoving 
volume is 



dn 
dM 



12 p dv 
TT M dM ' 



(3) 



where p is the present mean matter density. The peak 
height is defined as v{M) — 6c/cr{M), where a{M) is the 
linear rms fluctuation in a sphere containing a mean mass 
M . We compute cr{M) for an arbitrary cosmological model 
by integrating the linear power spectrum P\i-a{k) derived 
from the BBKS transfer functio n (Ba rdeen et al. 1986; 



with the conventions of Peacock 1997 ), evolved with the 
linear growth factor D{z). The density threshold 5c de- 
pends weakly on cosmology (i.e. on ri,„ and S^a) and was 
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Table 1. Cosmological Models 
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Fig. 3. Selection function (j){T, z) for ACDM model de- 
rived from the simulations. This function is the proba- 
bility that a cluster with temperature T and redshift z is 
detected and classified as an extended object in the survey 
catalogue. 
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Fig. 4. Selection function (f){S, z) expressed in terms of 
the flux S in the [0.5-2] keV band for the ACDM model. 
Only fluxes corresponding to temperatures above our 2 
keV limit are displayed. The vertical line shows a flux 
limit of 10"^'' erg s~^ cm~^; it illustrates deviations from 
a strict flux limit when realistic observing conditions are 
taken into account (esp. source confusion). 



computed using the fitting formulae of Kitiyama & Suto 



(1996) 
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4.2. Temperature Function 

The X-ray temperature of a cluster at redshift z is taken 
to be the virial temperature which is given by (see e.g. 
Eke, Cole & Frenk |1996|) 



kT 



7.75 



{i + z)ni 



M 



(—) 

V0.59/ 



\178 



keV, (4) 



where the average virial overdensity Ac{z, flm, ^a) can be 
evaluated using the fitting formulae of Kitiyama & Suto 
( |1996| ), Mi5 = IO'^^H-'^Mq and the value = 0.59 for 
the number of particles per proton mass corresponds to 
a hydrogen mass fraction of 76%. The factor f3y is equal 
to about 1 for a truncated singular isothermal sphere. We 
adopt this value as it provides a good fit to numerical 
simulations (Eke, Cole & Frenk 1996; Bryan & Norman 



1997) 



Combining Eqs. (||) and (Q) we can derive the differ- 
ential temperature function 



dn 



dn dM 
dM~dT' 



(5) 



It is often more convenient to consider the number den- 
sity of clusters with temperatures above a given minimum, 
ni>T)^J^dT"-lffl. 

To illustrate the dependence of our prediction on cos- 
mological parameters, we consider the three cosmologi- 
cal models hsted in Tab. i.e a tilted (r) CDM, ACDM 
and OCDM . The normalization of these models is de- 
termined by (Tsi the amplitude of mass fluctuations on 
8h~^ Mpc scale. Our chosen numerical values correspond 
to the constraints derived from current cluster surveys, 
namely as ~ 0.52r2~''-^^+''-^^^'" for the flat case and 
CT8 - 0.5 2a~°-^^+°-^°"'" for the for the open case (Eke et 



al. 199£). The shape of the matter power spectrum is con- 
trolled by the shape parameter F which, unless otherwise 
specified, we fix at 0.23, as i ndicat ed by galaxy clustering 
surveys (see Viana & Liddle 1996 and reference therein). 



4.3. Projected Cluster Counts 

From the temperature function ^ (Eq. ||^), wc can com- 
pute the projected surface density of clusters on the 
sky. Noting that the comoving volume element is dV = 
B?dxdVL^ where dSl is the infinitesimal solid angle, x is 
the comoving radius and R{x) is the comoving angular- 
diameter radius, we find that the number of clusters per 
unit solid angle, temperature and redshift interval is 



dN 
dTdz 



R 



,2 dx dn 
dz dT ' 



(6) 
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where g = --^ [(1 - n)a-'^ + fi^a-^ + Q^] as de- 
rived from the Friedmann equations. As a result the ob- 
served surface density of clusters per unit redshift 

dNobs f dN 



dT- 



-HT,z) 



(7) 



dz J dTdz ' 
where (j){T,z) is the survey selection function (Eq. |^). 

4.4. Predictions 

The predicted projected counts as a function of redshift 
are shown on Fig. |5| for the three cosmological models 
whose parameters are listed in Tab. |^. The counts cor- 
respond to the fuU 64 deg^ of the completed XMM-LSS 
survey. The counts from the three models agree at low 
redshifts (z < 0.2), as expected since the three models 
were normalized with the number of clusters in the lo- 
cal universe. On the other hand, the number counts dif- 
fer greatly at larger redshifts. The tCDM model predicts 
much smaller number of clusters at z > 0.2, while the 
ACDM and OCDM models differ for z > 0.6. The pre- 
dicted counts at z > 1 are larger in the OCDM model 
compared to that in the ACDM case, due to the somewhat 
slower evolution of the growth factor in the open model. 
These differences illustrate the well known fact that the 
evolution of cluster counts is a powerful probe of flm and 
Note that our predictions implicitly rely on the L — T 
relation of Eq. which was used to derive the selection 
function (p{T,z) from the simulations (see 

In the 64 deg^ of the survey, the expected number of 
detected clusters with < z < 1 (and T > 2 keV) is about 
900, 1000 and 175 for the ACDM , OCDM and tCDM 
model, respectively. For 1 < z < 2, the expected counts 
are about 400, 900 and < 10, for each model respectively. 
Our total number of clusters for < z < 2 of about 1300 
in the former model, is in approximate agreement with the 
predictions of Moscardini ct al. (2000) who predict about 
1100 clusters in the XMM-LSS for a ACDM model with 
slightly different flux and temperature thresholds. 

Fig. 1^ shows how these predictions for the ACDM 
model depend on erg and flm- Clearly, the number counts 
are very sensitive on these two parameters. Taking erg — 
0.93 ± 0.07 for fi^ = 0.3 (Eke et al. [l996[ ), the expected 
number of clusters in this model is between 600-1200 for 
< z < 1 and 200-700 for 1 < z < 2. 

The effect of the selection function on these predictions 
are shown in Fig. ^. The number counts for the ACDM 
model are shown as a function of redshift, for our selection 
function 4'{T, z) (as in Fig. ||), for the flux limited case and 
for the temperature limited case. For z ^ 0.9, the counts 
with the selection function are larger than the flux limited 
counts. This is expected, since a large number of clusters 
fainter than the flux limit contribute to the 4>{T, z) counts 
because of the tail of the selection function (see Fig. 
This demonstrates the importance of considering all the 
observational details when making such predictions, es- 
pecially confusion by point-like sources, an unavoidable 
drawback for a highly sensitive instrument. Notice that 





120 - 




100 - 


w 




00 




QJ 
X! 


80 '- 


^ 




CD 






60 - 






a, 




•z. 


40 - 




20 - 




0.0 



0.5 



1.0 

z 



1.5 



2.0 



Fig. 5. Projected number counts of clusters as a function 
of redshift in the three cosmological models. The selection 
function 0(T, z) for the XMM-LSS derived from image 
simulations was used for each model. The redshift bins 
have a width of Az — 0.05. 



the temperature limited counts agree with the other two 
for z < 0.4. This shows that our counts are limited mainly 
by the temperature at low redshifts. Not surprisingly, the 
T-limited counts overpredict the abundances for z > 0.4. 
As a check, we have compared our predictions with that 
of Romer et al. ( ^001 ). For their effective flux limit of 
'S'[o.5-2] > 1-5 X lO^i'^ erg s'^ cm'^ and kT > 2 keV, we 
find that our predictions agree very well with the counts 
shown in their Figure 6c. 

Note that the fluxes of the clusters measured in XMM- 
LSS will be rather uncertain (see VPG). This will prevent 
us from making accurate determinations of the luminos- 
ity, and therefore of the temperature and of the mass, of 
each cluster. In this paper, we therefore consider statis- 
tical quantities (such as the counts and the correlation 
function) which are averaged over the population of clus- 
ters detectable in the survey, and therefore do not require 
this information. A study of the further constraints which 
can be derived from the (uncertain) flux measurements is 
left to future work. 



4.5. Cosmological Constraints from Cluster Counts 

We now wish to study how the cluster counts can be used 
to constrain cosmological parameters. For this purpose, 
we generated cluster counts from the predicted counts 
from Eq. (^ in ACDM model for several redshift bins 
(as shown in Fig. ^ for different models). We used the 
XLSS selection function 0(r, z) and survey area, and con- 
sidered redshift bins of width Az = 0.1 from z = to 
z = 1 assuming Poissonian errors for the number counts 
in each bin. Note that the Poisson errors neglect the uncer- 
tainty arising from cosmic variance; this is a good approx- 
imation since this latter uncertainty is small for a survey 
area as large as that of the XLSS (compare the correla- 



Refregier et al: Clusters in the XMM-LSS 

10000 



100 



80 



T > 2 keV 
Si(T,z) 

S. =10"" erg/s/cm^ 




2.0 



Fig. 6. Projected number counts of clusters for the ACDM 
model with different selection criteria: with the selection 
function 4>{T, z) (as in figure with a temperature limit 
only (fcT > 2 keV), and with a flux limit (5[o.5_2] > 10"^'' 
^). Again, the redshift bins have a width of 



erg s" 
Az = 



cm 
0.05. 



tion length ^ 20h~^ Mpc to the survey size, ^ 200h~^ 
Mpc at z ~ 0.5). We then calculated the x^-difference be- 
tween the simulated counts and that expected in a ACDM 
model for a range of values for o-g and flm ■ This approach 
is rather computer intensive, but it is more accurate than 
the faster Fisher matrix method, which approximates the 
confidence regions as ellipses and can be incorrect when 
the parameter space is non-Gaussian (Holder et al. 2001 ). 
Note that our fitting method is only strictly valid for 
gaussian errors. This is however a good approximation in 
our case, given the large number of clusters in the redshift 
bins shown in Figure |^, further enhanced by our use of 
Az = 0.1 bins rather than 0.05 in the Figure. 

The resulting constraints on these cosmological pa- 
rameters are shown as the solid lines in Fig. |^. The in- 
put ACDM model is that of Tab. |l| and is shown with 
a cross. The 90% and 95% confidence regions are shown 
as solid lines. The cluster counts alone will provide tight 
constraints on both parameters, with 95% uncertainties 
of about 0.06 and 0.05 for erg and ^m, respectively. As 
can be seen from the elongation of the contours, the two 
parameters are however somewhat degenerate. Any addi- 
tional information on either parameters, can thus be used 
to reduce their respective uncertainty. Note that these 
constraints only reflect statistical errors, and neglect po- 
tentially important systematic uncertainties in the scaling 
laws of clusters. A discussion of the limitations imposed 
by systematics on our predictions is presented in §|^. 

It is also interesting to investigate how cluster counts 
can constrain F, the shape parameter of the matter power 
spectrum. For this purpose, we follow the same Pro- 
cedure, this time varying F and flm, while erg follows the 
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Fig. 7. Cluster counts expected for the XMM-LSS as a 
function of CTg and SI™ in the ACDM model. The XMM- 
LSS simulated selection function was used in all cases. 
Counts for the < z < 1 and 1 < z < 2 redshift intervals 
are shown as the solid and dashed lines, respectively. In 
each case, models with = 0.4, 0.3 and 0.2, are shown 
from top to bottom, respectively. 



Eke et al. (|l996| ) relation (see §4^). The resulting confi- 
dence contours for the ACDM model are shown as the solid 
lines in Fig. ||. Clearly, these two parameters are quite de- 
generate with counts alone, hampering the determination 
of r2„i, which can only be determined with an accuracy 
of about 40% (95% CL). More information is therefore 
required to alleviate these limitations. One obvious pos- 
sibility is to use other measures of large-scale structure 
such as galaxy catalogues to constrain F. This has the 
disadvantage of relying on assumptions about the bias of 
galaxies and on an external data set. In the next section, 
we show how the degeneracy can be broken by measuring 
the correlation function of the galaxy cluster population. 



5. Cluster Correlation function 



We now turn to the cluster correlation function which 
quantifies the spatial clustering of clusters in the survey. 
We first use the extended halo model of Mo & White 
( 1996| ) to predict the observed cluster correlation function 
at a given redshift. We then consider the average correla- 
tion function in a finite redshift range, compute its uncer- 
tainties, and study the constraints its measurement can 
place on cosmological parameters. 
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Fig. 8. Constraints on the cosmological parameters ftm 
and (Tg for the ACDM model, from the cluster counts (solid 
lines) and from the correlation function (dashed lines) . In 
each case, the 68%, 90% and 95% confidence level contours 
are shown along with the assumed model (cross). 
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Fig. 9. Constraints on F and fim for the ACDM , from 
the cluster counts (solid lines) and from the correlation 
function (dashed lines). As before, the 68%, 90% and 95% 
confidence level contours are shown (solid lines) along with 
the assumed model (cross). 

5.1. Observed Correlation Function : at a single 
redshift 



According to the Mo & White (1996) halo formalism, the 
correlation function of two sets of clusters with masses M 
and M' and with comoving separation r is given by 

C(r, M, M', z) ~ h{M, z)b{M', z)^Ur, z), (8) 



where ^iin(r, z) = ^iin(r, O)Z)^(z) is the linear correlation 
function which is the Fourier Transform of the linear 
power spectrum Piin(fc,z). The mass dependent bias pa- 
rameter of the halos is given b{M, z) = 1 + (v^ — l)/(5c, 
with the conventions of §4.1. Note that the separation of 



the clusters is assumed to be small compared to the scale 
in which any evolution takes place. 

It is easy to show that the resulting observed correla- 
tion function in a narrow redshift interval is given by (Suto 
et al. 200C; Moscardini et al. 200C and reference therein) 



fobs(?-,2) = 6^ff(z)^iin(r,z), 



where the effective bias is 



dM^^b{M) 



dM 



dM- 



dn, 



obs 



dM 



(9) 



(10) 



and where the observed differential number counts are 
given by ^^^^ = ^j^4>iT, z). Note that, in our analysis, 
we neglect redshift-space distortions which were shown to 
yield only about 10% corrections on th e amplitude of the 
corre lation function (Suto et al. 2000 ; Moscardini et al. 
^OOOD 

In general, the evolution of ^obs(?', z) is determined by 
two competing effects. First, the growth of structures in- 
duces the linear correlation function ^ii„(r, z) to decrease 
as the redshift increases. On the other hand, the clus- 
ters which are detectable at large redshifts are more mas- 
sive and therefore more strongly biased. The effective bias 
&cff(z) thus tends to be larger at high redshift. To study 
the interplay between these two effects it is convenient to 
define 



bcS,o{z) 



Cobs(r, z) 



6in(r,0) 



b,s{z)D{z). 



(11) 



This quantity provides the bias of the observed cluster 
correlation with respect to the linear correlation function 
at z = 0, and therefore quantifies the evolution of the 
correlation function. 

The behaviour of beSfi{z) is shown on Fig. |l^ for the 
three cosmological models and for the XMM-LSS selection 
function. The curves are remarkably flat for all models, 
showing that the cluster correlation function evolves only 
very weakly from z = to z = 2. It is interesting to study 
whether this lack of evolution depends on the selection 
function. Fig. |l^ shows 6eff,o(-2^) for each selection scheme 
for the ACDM model. The evolution is also very weak for 
the temperature limited case. For the flux-limited sample, 
the evolution is stronger for z > 1. This evolution is lost 
when the full selection function 4>{T, z) is used. In all cases, 
there is effectively no evolution for < z < 1. We shall 
thus, in the following, derive the constraints on cosmology 
integrating the correlation function over the < z < 1 
range. 

5.2. Observed Correlation Function: Redshift Average 

To maximize the sensitivity, it is useful to measure the cor- 
relation function over a wide redshift range. In this case. 
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Fig. 10. Evolution of 6cff,o(z), the effective bias of the 
cluster correlation function ^obs('', 2) with respect to the 
linear mass correlation function ^iin(?', 0) at z = 0. This 
quantifies the evolution of the cluster correlation function. 
It is shown for the three cosmological models as a function 
of redshift. 
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Fig. 12. Prediction for the cluster correlation function av- 
eraged in redshift interval < z < 1. The three cosmolog- 
ical models are shown. For clarity, the expected measure- 
ment error bars for the full 64 deg^ of the XMM-LSS are 
shown only for the ACDM model. . 
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Fig. 11. Evolution of 6off,o(-z) in the ACDM model for dif- 
ferent selection schemes: the XMM-LSS selection function, 
a temperature limited survey, and a flux limited survey. 



we must take into account the pair-weighted evolution of 
the correlation function ^ (r, z) within the light-cone sec- 
tion of interest. The average correlation function in the 
redshift interval Zmin < z < z„iax is thus given by (Suto et 
al. 2000| ; Moscardini et al. pOOO and reference therein) 



(12) 



where the observed number density of clusters is nobs('Z) 
/ dTj^<l>{T,z). 



The resulting correlation functions for the three cos- 
mological models are shown on Fig. ^ The XMM-LSS 
selection function was used in all cases, along with a red- 
shift range of < z < 1. The correlation functions for the 
3 models have very similar shape, as expected since the 
same value for F was assumed in all cases. Notice that the 
tCDM correlation function has a lower amplitude than 
that for the other two models. This is expected since the 
former model has a lower value for as (see Tab. The 
ACDM and OCDM models have very similar values of 
cTg and thus yield correlation functions with very similar 
amplitudes. 

Our predictions a re in g ood agreement with the results 
of Moscardini et al. 



2000 



In particular, for our ACDM 
model, these authors find a correlation length tq (defined 
by lobsiro) = 1) of about 13 and 15 h^^ Mpc for z < 0.3 
and z > 0.3, respectively. This is to be compared with 
our value of tq ~ 17 h~^ Mpc for 2: < 1. We also verified 
that the evolution of the correlation function is weak if we 
adopt their assumed strict flux limit of S'[o.5-2] = 5 x 10^^^ 
erg s~^ cm~^ . 



5.3. Uncertainties 

We now wish to estimate the uncertainty involved in mea- 
suring the cluster correlation function. The simplest way 
to measure the correlation function is to count the number 
of pairs in excess of random in the survey volume. (This is 
sometimes taken as the definition of the correlation func- 
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tion). An estimator for $obs('') ^ bin of radii between r 
and r + Ar is thus 



Usir^Ar] 



jyohs 



pairs 



(r, Ar) 



pairs 



(r, Ar) 



(13) 



where -/Vpafrs and -/Vpafrs are the number of pairs (i.e. with 
clustering) and for a random distribution (i.e. without 
clustering), respectively, in the observed radius interval. 
Ignoring boundary effects, the random number of pairs is 
given by 



(14) 



where A is the solid angle of the survey. 

For weak signals, the error in measuring ^^^^g(r,Ar) 
will be dominated by Poi sson s tatistics, and will thus be 
given by (see e.g. Peebles 1980 ) 



(15) 



This provides us with an expression for the error in the 
correlation function for a finite survey. The resulting errors 
for the ACDM model and for the full 64 deg^ of the XMM- 
LSS are shown in Fig. 12. Again, these errors only reflect 
statistical errors and neglect uncertainties in the cluster 
scaling relations. A discussion of the impact of systemat- 
ics, including that arising from spatial variations of the 
survey sensitivity, is presented in §|[ 

5.4. Combined Cosmological Constraints 

We now study how the measurement of the cluster corre- 
lation function constrains cosmological parameters. As for 
the cluster counts (§ |3.4[ ), we use a x^-fit to simulated mea- 
surements of the correlation function to derive confidence 
regions in parameter space. We considered a redshift inter- 
val of < z < 1 and computed the errors using Eq. (|l5|). 

The resulting confidence contours for a joint fit of 
and (78 are shown as the dashed lines in Fig. ||, for the 
ACDM model. The constraints from the correlation func- 
tion are rather weak on this plane alone, with little de- 
pendence on erg ■ These constraints are however somewhat 
orthogonal to that from cluster counts and are thus com- 
plementary. 

The constraints from the correlation function on F and 
are shown in Figure ^. The shape of the resulting confi- 
dence contours agrees qualitatively with the predictions of 
Moscardini et al. ( 2000) for the shallower REFLEX survey. 
The contours show little dependence on ftm and are there- 
fore nearly orthogonal to that from cluster counts alone. 
With the combined counts and correlation function con- 
straints, and F can be measured with a precision of 
about 15%, 10%, and 35%, respectively at the 95% confi- 
dence level. The cluster correlation function thus provides 
the required information to break the degeneracy present 
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Fig. 13. Constraints on flm and as resulting from a fit 
of the OCDM model to the ACDM predictions. Both the 
cluster counts (solid lines) and the correlation function 
(dashed lines) are shown. In each case, the 68%, 90% and 
95% confidence level contours are shown, along with the 
minimum (cross) and the input values for the ACDM 
model (filled triangle). 



when cluster counts alone are considered. The impact of 
systematics on these constraints is discussed in §|[ 

Until now, we have focused on the ACDM model, thus 
implicitely assuming prior knowledge that the universe is 
fiat. While flatness is strongly indicated by measurements 
of the Cosmic Microwave Background anisotropics (e.g. 
Jaffc et al. [2001 ), it is interesting to establish whether our 
predictions depend on this assumption. To study this, we 
fitted an OCDM model to the ACDM predictions for both 
the cluster counts and the correlation function. The result- 
ing constraints on as and Qm are shown in Fig. |l^. For the 
cluster counts, the best fit values for both of these param- 
eters are now biased (compare the best fit, cross, to the 
input values, triangle). This bias again hampers the de- 
termination of using cluster counts alone. Thankfully, 
the correlation function can again help overcome this limi- 
tation. Indeed, the contours for i^(r) are now somewhat in- 
compatible with that from the counts (compare to Fig. ||) . 
The consistency between the counts and correlation func- 
tion constraints can thus be used as a diagnostic and as a 
discriminant between the ACDM and the OCDM model. 



6. Systematic Uncertainties 

So far in our analysis, we have neglected systematic uncer- 
tainties. Here we review the major source of systematics 
and discuss their impact on our predictions.^ 



^ We thank Kathy Romer, the referee, for her many sugges- 
tions which are included in this section. 
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First, we have assumed that the L ~ T and M — T 
relations were perfectly known. Physical processes such 
as feedback from star formation and reionization, can 
modify these relationships and therefore yield detectable 
changes in the cluster counts and correlation function (see 
e.g. Muanwong et al. 2001 and reference therein). These 
could perhaps explain the discrepancy which has recently 
been found between the observed M — T relation and 
that derived from numerical simulations (see Finoguenov, 
Reiprich, & Bohringer 2001; Allen, Schmidt, & Fabian 
2001 and reference therein). Cluster mergers have also 
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been shown to produce sharp temperature and luminosity 
jumps, thus likely affecting the L — T relation (Ricker & 
Sarazin 2001 ). Another important issue is the existence of 
high redshift cooling flows which may affect the detectabil- 
ity of distant clusters and their L — T relation (e.g. Henry 
2000| ). 

By the time the XMM-LSS is completed, better in- 
sights will be gained about these various physical effects, 
thanks to deep pointed observations of clusters with XMM 
and Chandra. An interesting way of parameterizing the 
residual uncertainties has been proposed by Diego et al. 
(2001). These authors kept the normalisation and evolu- 
tion of the different cluster scaling relations as free pa- 
rameters, and studied the resulting degeneracies with cos- 
mological parameters. It would be interesting to extend 
their approach to include constraints from the cluster cor- 
relation function. The study of these pending theoretical 
questions is left for future work. 

In our simulation, we have also assumed that the core 
radius Vc is fixed. More realistically, the core radius could 
vary w ith mass and redshift. For instance, Jones et al. 
(|1999D found a correlation between core radii and temper- 
atures in a cluster sample from the Einstein mission. Our 
detection sensitivity depends most strongly on the clus- 
ter flux, which is independent of the core radius in our 
simulations. A varying Tc is therefore unlikely to affect 
our selection function very much, except perhaps in our 
ability to separate clusters from point sources. 

Another systematic effect arises from our assumption 
of a sharp temperature cutoff [kT > 2 keV) for our cluster 
selection. Since temperatures will not be available with 
XMM-LSS, this will be enforced in practice by setting a 
luminosity cutoff, derived from the observed fluxes and 
redshifts. Because the luminosity function is very steep, 
the errors in the observed fluxes will lead to a bias tending 
to include low mass clusters in the counts. This effect can 
be corrected for by running simulations similar to ours, 
but which would include objects with kT < 2 keV. 

Another complication is induced by variations in the 
detection sensitivity over the survey area. Figure |l^ shows 
the effective exposure time for the XMM-LSS tiling strat- 
egy, in which the 10 ksec exposures are separated by 20 
arcmin. As is apparent in the figure, vignetting causes 
variations in the exposure time of about 50%, correspond- 
ing to flux sensitivity variations of about 25%. This will 
affect both the uniformity of the cluster counts and the 
observed cluster correlation function. Since the vignetting 




Fig. 14. Exposure time map for the XMM-LSS tiling pat- 
tern. Four XMM fields are displayed out to a maximum 
radius of 13'. The fields are separated by 20 arcmins from 
each other. 

function is well calibrated, the sensitivity variations can 
however be accurately predicted and corrected for. This 
can be done by using further image simulations to com- 
pute the selection function as a function of position. Note 
that, in the simulations we have presented in we have 
only computed the average selection function by placing 
clusters throughout the inner 13' of a single XMM-EPIC 
pointing. Our detection limits are thus conservative, as 
they do not include the increased sensitivity afforded in 
the regions of overlap of two XMM fields. 

7. Conclusions 



Following the R EFLE X survey (Guzzo et al. |1999 , 
Bohringer et al. ^OOl] ) based on the ROSAT All-Sky- 
Survey and the NEP ROSAT survey (Henry et al. ^001]) , 
the XMM-LSS survey will be about 1000 and 10 times 
deeper, respectively, and thus opens wide prospects for 
cosmology. Indeed, it will provide an independent mea- 
surement of cosmological parameter and thus comple- 
ments Cosmic Microwave Background and Supernova ex- 
periments. It will also provide a test of important ingre- 
dients of the standard cosmological model, such as the 
gravitational instability paradigm and the gaussianity of 
initial fluctuations. 

To study the constraints XMM-LSS will set on cosmo- 
logical parameters, we first derived the selection function 
for the survey using detailed simulations of cluster detec- 
tion in XMM-Newton images. We found that our selection 
function differs significantly from a simple flux-limit selec- 
tion. We then computed the expected number counts of 
clusters in several CDM models. We found that, for the 
currently favoured ACDM model, we expect about 600- 
1200 clusters to be detectable in XMM-LSS at < z < 1 
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and about 200-700 at 1 < z < 2, the uncertainty being 
dominated by the current errors on erg and f2m. 

Clusters counts beyond z > 0.5 depend on cosmologi- 
cal parameters and can thus be used to constrain models. 
Within a ACDM model, we found that the rcdshift de- 
pendence of the XMM-LSS cluster counts will allow us 
to measure erg and with a precision of about 6% and 
18% (95% CL), respectively, if the shape parameter T is 
known. In the absence of prior knowledge on F, the pre- 
cision on these parameters degrades considerably, if only 
cluster counts are considered. 

This limitation can be circumvented by considering 
the cluster correlation function. One of the strength of 
XMM-LSS is indeed to offer a uniform coverage over a 
wide contiguous area (64 deg^), with an extensive spectro- 
scopic follow-up. This will allow us to measure the corre- 
lation function in several rcdshift bins out to z = 1. Using 
the extended halo formalism of Mo & White ( 1996 ), we 
computed the correlation function of clusters detectable 
in XMM-LSS. We find that, for the selection function of 
XMM-LSS, the correlation function is not expected to 
evolve significantly from z = to 2. This results from 
the competing effects of the growth of mass perturbations 
and the stronger bias of the detectable massive clusters 
at large redshifts. It will thus be important to compare 
the correlation function measured in two rcdshift bins be- 
tween z — and z = 1: the verification of the lack of 
evolution provides a test of the bias model for haloes and 
of the gravitational instability paradigm. 

The amplitude and shape of the cluster correlation 
function can be used to lift the degeneracies present when 
cluster counts alone are considered. Within a ACDM 
model, the correlation function function measured in 
XMM-LSS at < z < 1, combined with the cluster counts, 
will constrain fl„n erg and T with a precision of about 15%, 
10% and 35%, respectively (95% CL). Moreover, the com- 
bination of the counts and the correlation function will 
provide a consistency check for the ACDM model, and a 
discrimation between this model and the OCDM model. 
The XMM-LSS therefore has great potentials for the mea- 
surement of cosmological parameters. Note however that 
the above constraints only reflect statistical errors and ig- 
nore the limitations arising from systematic uncertainties. 
Instrumental systematics, such as sensitivity variations, 
could be important but can be accurately corrected us- 
ing further image simulations. Systematic uncertainties in 
the cluster scaling relations may also contribute signifi- 
cantly to the error budget. Bv the time the XMM-LSS 
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